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ABSTRACT 


This  thesis  investigates  the  synthesis  of  smoothed  impulse  responses  using  sampled 
data  of  the  input  and  output  of  random  noise  driven  electromagnetic  systems.  Special 
interactive  software  was  developed  for  NPS's  time  domain  electromagnetic  scattering 
laboratory.  The  system  performs  signal  acquisition,  synthesizes  time  and  frequency  do¬ 
main  scattering  responses  using  broad  band  random  noise  and  provides  results  as  easily 
evaluated  graphic  displays.  Attempted  validations  of  the  system  are  made  by  computing 
synthesized  impulse  responses  for  microwave  filters  and  transient  scattercrs  to  alternate 
experimental  and  computational  data. 
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I.  INTRODUCTION 


A.  OVERVIEW 

The  objective  of  this  research  is  to  demonstrate  the  viability  of  performing  high- 
resolution  impulse  response  scattering  measurements  using  a  broad-band  noise  source. 
The  development  of  laboratory  facilities  for  high-resolution  impulse  response  scattering 
measurements  has  generally  proceeded  using  two  major  techniques.  The  first  of  these 
employs  a  stepped-frequency  coherent  oscillator  and  a  vector  (magnitude  and  phase) 
receiver.  Impulse  response  target  characteristics  are  obtained  via  inverse  Fourier  trans¬ 
formation  of  the  frequency  domain  data.  The  second  method  obtains  impulse  response 
measurements  directly  in  the  time-domain  by  use  of  a  repetitive  fast-pulse  target  illumi¬ 
nation  with  a  sampling-scope  acting  as  the  receiver.  The  time-domain  approach  for 
scattering  measurements  offers  a  viable  alternative  to  the  more  prevalent  continuous 
wave  approach.  Transient  scattering  measurements  provide  waveforms  that  can  be  more 
directly  interpreted  as  to  cause  and  effect  and  allow  exact  range-gating  of  target  re¬ 
sponses  for  elimination  of  unwanted  clutter  (lief.  1.  2]. 

The  development  of  digital  sampling  oscilloscopes  and  broad  band  noise  generators 
made  possible  a  third  method  for  high-resolution  impulse  response  scattering  measure¬ 
ments.  This  technique  is  based  on  the  estimation  of  the  crosscorrelation  function  be¬ 
tween  input  and  output  using  a  broad-band  noise  source  as  the  transmitter.  Ibis 
methodology  has  been  employed  in  measuring  the  mechanical  impulse  response  of  large 
structures,  such  as  bridges  and  buildings  [Ref.  3). 

The  practical  advantages  of  this  third  technique  are  two-fold.  The  first  is  the  wide 
bandwidth  and  high  power  available  from  noise  sources  which  are  available  at  much 
lower  cost  than  comparable  stepped-frequency  sources.  A  conclusive  demonstration  of 
the  viability  of  noise  source  impulse  response  measurements  may  lead  to  further  refine¬ 
ments  and  ultimately  to  commercial  marketing  of  this  technology.  The  second  advan¬ 
tage  is  related  to  the  use  of  noise-source  illumination  for  tactical  and  strategic  radar 
applications.  An  obvious  benefit  would  be  the  masking  of  the  radar  interrogation  signal; 
this  would  appear  at  the  target  as  either  jamming  or  interference.  Confusion  would  re¬ 
sult  as  to  proper  countermeasures  to  be  employed  by  the  target.  Fuitlicrmore.  since  this 
method  uses  random  noise,  which  is  orthogonal  to  other  signals  and  other  noise 
(producing  zero  crosscorrelation  with  these),  a  high  SNR  requirement  could  be  achieved. 


This  thesis  is  a  continuation  in  a  series  of  efforts  in  transient  electromagnetic  scat¬ 
tering  that  began  in  1979  at  the  Naval  Postgraduate  School  (NPS).  The  development 
of  digital  sampling  techniques  allowed  the  development  of  transient  scattering  ranges  in 
the  late  1960's  [Ref.  4],  A  transient  scattering  range  having  sufficient  bandwidth  and 
signal  to  noise  ratio  (SNR)  to  support  radar  target  identification  was  initially  con¬ 
structed  at  NPS  in  1980  using  a  ground  plane  configuration  [Ref.  5,6].  A  very  wideband 
(short  pulse)  free-field  scattering  range  was  then  constructed  in  1983  for  implementation 
of  target  identification  based  upon  natural  resonances,  as  introduced  by  Mittra  and  Van 
lilaricum  [Ref.  7|.  This  free  field  transient  range  was  validated  by  Mariategui.  and 
McDaniel  [Ref.  8.9]. 

3.  NOISE  SOURCE  IMPULSE  RESPONSE  MEASUREMENT  CONCEPT 

In  the  noise  source  impulse  response  measurement  technique,  a  broad  band  noise 
signal  is  used  as  the  input.  The  system  impulse  response  is  derived  from  measuring  the 
crosscorrelation  function  between  the  input  and  output. 

The  configuration  shown  in  Figure  la  has  been  used  to  implement  the  noise  source 
impulse  response  measurement  technique  in  the  NT’S  scattering  laboratory.  It  was 
thought  at  first  that  the  Tektronix  7804  Digital  Processing  Oscilloscopes  (DI’O)  would 
allow  simultaneous  measurements  of  both  channels.  After  much  effort,  this  was  found 
to  be  untrue,  as  will  be  shown.  A  new  method  was  developed,  using  the  measurement 
of  the  difference  of  the  input  and  output  time  sequence.  A  digital  estimator  of  the 
crosscorrelation  function  is  formed  in  the  data  processor  for  the  measurement  and  the 
expectation  of  the  crosscorrelation  function.  R„(n),  is  obtained  by  further  averaging  of 
this  estimator  using  N-block  averages  of  the  stored  data. 

The  complete  scattering  measurement  using  a  noise-source  requires  three  steps. 
Looking  at  the  scattering  range  system  representation  in  Figure  lb.  the  target  transfer 
function  //,(/)  will  be  evaluated  by  measuring  X(n)  and  I'(«).  the  input  and  output  of  the 
respective  transmitting  and  receiving  antennas.  The  second  step  necessitates  factoring 
out  the  effects  of  the  antennas,  in  addition  to  elmiinating  the  signal  pollution  of  the 
spurious  cross-coupling  and  clutter  in  the  chamber.  The  following  three  measurements 
ate  made  with  equivalent  system  diagrams  shown  in  Figure  2: 


Step  1.  :  Background  measurement  (no  target) 


clutter 

__ , 

L_LJ 

J 

i -  [_ 

Trans.  Ant. 

Hc  - 

Rec.  Ant. 

Step  2. 


Hi(f)  =  Ht(f)  Hc(f)  Hr(f) 

Calibration  target  (sphere)  measurement 


Step  3. 


H2(f)  =  Hi(f)  +  Ht(f)  Hp(f)  Hr(f) 
Test  target  measurement 


H3(f)  =  Hi(f)  +  Ht(f)  H»(f)  Hr(f) 

Step  4.  :  Subtraction  of  the  background  effect 
H«(f)  =  H2(f)  -  Hi(f) 

Ha(f)  =  Ha(f)  -  Hi(f) 

Step  5.  :  Extraction  of  the  target  response  (Deconvolution) 

Hs(f)  _  Ht(f)  Ha(f)  Hr(f)  _  H»(f) 

H«(f)  Ht(f)  Hp(f)  Hr(f)  Hp(f) 

For  optimal  compensation,  use  Riad's  method 

Ha(f) _ Ha(f)  H<(f) _ 

H4(f)  ~  |  H-t(f)  |  *  +  (£  H«*(f))/N 

Therefore, 

H  (f\  —  Hs(f)  H«(f)  .. 

~  |  H<(f)  J 2  +  (S  H«2(f))/N  Hpc(f) 


Hpc(f) 


Hpc(f)  :  computed  frequence  response  of  the 
calibration  target,  (sphere) 


Figure  2.  The  Procedure  of  Noise  Source  Impulse  Response  Scattering  Measure¬ 
ment 


1.  No  target  present,  measure  /i,(r) 

2.  Calibration  target  (sphere)  present,  measure  /i:(r) 

3.  Test  target  present,  measure  /i,(t) 

Then,  A,(r)  will  be  subtracted  from  h2(t)  and  /«,(/)  to  eliminate  the  clutter  e fleet. 

/i4(/)  =  /t2(r)-/j,(r)  (1.1) 

/':(')  -/»j(/) -/»((/)  I  1.2) 

Finally,  the  desired  impulse  response  is  extracted  from  /i4( /)  and  fu(t)  using  an  optimal 
deconvolution  technique  known  as  Riad's  method  [Ref.  R»|  which  will  be  described  in 
detail  in  the  following  chapter. 

An  important  difficulty  encountered  in  this  thesis  research  was  that  no  sampling 
device  with  a  high  enough  rate  for  the  source  signal,  having  a  bandwidth  of  approxi¬ 
mately  13  GHz.  is  available  thus  far.  Consequently,  the  estimation  of  the  crosscorrc- 
laiion  function  could  not  be  computed  using  a  pair  of  propet  ly  sampled  sets  of  the  time 
function  of  each  input  and  output  signal.  Two  alternate  methods  were  developed  which, 
for  the  case  of  ergodic  random  noise,  permits  arbitrarily  slow  sampling  to  be  used.  These 
two  Nyquist  independent  techniques  will  be  described  in  Chapter  II. 

C.  OVERV  IEW  OF  THESIS 

The  objectives  of  this  thesis  were  to: 

1.  Develop  a  working  software  program  which  will  acquire  transient  response  data 
from  the  target  and  compute  the  estimation  of  the  crosscorrelation  function  fol¬ 
lowed  by  a  computation  of  the  system  impulse  response  using  a  deconvolution 
technique. 

2.  Demonstrate  impulse  scattering  response  measurements  using  a  broad  band  ran¬ 
dom  noise  source,  and  verify  the  performance  in  comparison  with  another  tech¬ 
nique:  time-domain  measurement  using  modified  step  function  input. 

Chapter  II  will  expand  the  theory  of  noise  source  impulse  response  measurement. 
This  will  examine  the  basic  theory  of  analog  crosscot  relation  techniques  for  acquiring 
system  impulse  responses  and  consider  an  appropriate  derivation  for  the  discrete  version 
of  the  technique  for  the  sampled  signal.  In  addition,  it  will  quantify  the  estimation  of 
the  crosscorrelation  function. 

Chapter  III  describes  the  experimental  system  in  detail.  I  his  will  include  the  de¬ 
scription  of  the  original  Digital  Processing  Oscilloscope  (DPO)  and  modifications  that 


were  made.  A  description  of  the  noise  source  hardware  and  its  modified  hardware  will 
also  be  given  in  this  chapter.  In  addition,  the  problems  encountered  in  this  research  will 
be  discussed  and  some  guidelines  for  the  laboratory  work  will  be  summarized. 

Chapter  IV  contains  the  calibration  and  validation  measurements  for  simulated 
targets.  This  will  include  the  initial  quality  tests  of  simultaneous  channel  sampling,  noise 
source,  and  crosscorrelation  function  that  were  made  by  sampled  data  using  the  noise 
source.  The  derivation  of  the  crosscorrelation  of  the  measurement  itself 
(autocorrelation)  and  frequency  response  will  be  considered  using  TFT  techniques.  The 
inirial  test  involving  a  microwave  bandpass  amplifier  using  a  noise  source  and  a  step 
generator  with  amplifier  will  also  be  discussed. 

Chapter  V  will  describe  the  electromagnetic  scattering  measurement.  The  exper¬ 
imental  setup,  scattering  range  and  characteristics  of  the  antenna  used  will  be  explained. 
Finally,  the  impulse  response  measurement  of  a  metal  sphere  target  will  be  attempted 
using  the  noise  source.  The  resultant  failure  of  this  measurement  will  be  considered. 

Chapter  VI  provides  some  conclusions  about  this  experimental  appion.h  involving 
noise  source  impulse  response  measurements.  Additionally,  recommendations  arc  made 
regarding  improvements  in  the  current  system  and  possibilities  for  further  icscarch. 


II.  THEORY  OF  NOISE  SOURCE  IMPULSE  RESPONSE 

MEASUREMENT 


A.  THEORY  OF  ANALOG  SIGNAL  CROSSCORRELATION 

Looking  at  the  system  diagram  in  Figure  lb,  the  responses  of  transmitter  antenna 
h,(t )  and  receiver  antenna  h,(t)  are  cascaded  with  the  parallel  responses  of  the  target  h,(t) 
and  the  clutter  ftc(r).  forming  an  overall  response,  h{i), 

h(t)~  hXD  *  [  h.(t)  +  /ic(0  ]  *  /[,(/)  (2.1) 

The  total  measurement  system  can  be  represented  as  a  simple  linear  system  whose  re¬ 
sponse,  1(/),  due  to  an  input,  A'(0.  is  given  by  the  convolution 

F(0  =  /i(0  *  AT/) 

i "  ’) 

=  h{a).\(t  ~  a)  da 

Jo 

It  can  be  shown  that  the  crosscorrelation  of  the  input  and  output  of  a  linear  system 
estimates  the  system  impulse  response  when  the  input  has  a  bandwidth  that  is  large 
compared  to  the  bandwidth  of  the  system.  This  results  from  superposition  as  applied 
to  stochastic  expectation, 

Rxy(t)~E[X(r)nr  +  ,)\ 

=  £(/«(/)  *A'(r)A'(r  + 01 
=  /,(/)*£(  A'(r)A'(r  +  0) 


Rxy(t)  =  hit)  *  R,(0 


(2.3) 


For  crgodic  processes,  the  mean  values  and  moments  can  be  determined  by  time  aver¬ 
ages  as  well  as  by  ensemble  averages. 

RxyV)  =  <  A'(r)K(r  +  /)> 


=  r  Mt)Y(t  +  i)Jz 

J  — >5 


(2.4) 


Replacing  the  output  signal  K(/)  by  the  convolution  defined  in  Equation  (2.2)  results  in 


Figure  3.  Crosscorrelation  Measurement  System  Model 


This  leads  to  an  important  property  of  the  Fourier  transform  relation,  which  is 
SxyW  =  H{f)  Sx(f>  (2.6) 

known  as  the  crosscorrelation  theorem  [Ref.  1 1|. 

B.  CROSSCORRELATION  MEASUREMENT 

1.  Method  1. :  Simultaneous  Dual  Channel  Measurement  with  Sampling  Rate 
Greater  than  Nvquist  Frequency. 

Three  methods  of  performing  the  measurement  could  be  used  to  estimate  the 
crosscorrelation  of  the  linear  system.  One  possible,  probably  the  most  general,  way  is 
to  directly  sample  the  input  and  output  sequence  simultaneously  and  store  the  data  in 
the  computer  memory.  This  involves  a  delay  of  the  input  sequence,  and  a  computation 
of  the  average  of  the  product  of  the  two  sequence  vectors. 


This  scheme  is  illustated  in  Figure  4  and  compared  to  other  methods.  Two  conditions 
must  be  satisfied  to  use  this  method. 

1.  Sampling  must  be  done  with  a  sampling  rate  greater  than  the  N'yquist  frequency, 

2.  Sampling  of  two  channels  must  be  done  simultaneously. 

2.  Method  2. :  Simultaneous  Dual  Channel  Measurement  with  Arbitrary  Sampling 
Rate.  (Less  than  Nyquist  Frequency) 

Method  1  could  not  be  used  because  the  required  bandwidth  of  the  DPO 
sampler  is  insufficient.  An  alternative  method  requires  manually  shifting  the  input  (or 
output)  signal,  being  incremented  by  the  sampling  interval.  Samples  of  the  input  and 
ou!  put  signal  are  taken  with  a  suitable  sampling  rate,  since  each  sampled  time  sequence 
is  treated  as  the  sample  set  of  the  ensemble  space  rather  than  a  time  sequence. 

AT- 1 

Rxy(n)  ~  ^  i+nl'.k 

k= o 

A  — 1 

R.ry(")  =  ~7  y.\ XkL)  Y(kL  +  n)  (2.8) 

k=0 

where, 

/t  =  «T  with  ~  >  2/max  (2.‘>) 

T.  s  LT  ,  L$>  l  (2.10) 

I  lei c.  T  denotes  the  time  period  of  delay  and  Ts  is  the  sampling  period.  Using  this 
technique,  the  estimation  of  the  crosscorrclation  of  the  delayed  time  point  can  be  com¬ 
puted  by  each  measurement  of  input  and  output  signal  with  an  appropriate  amount  of 
delay.  This  method  is  illustrated  in  Figure  4.  l  ive  remaining  condition  which  must  be 
satisfied  for  the  crosscorrelation  measurement  is  simultaneous  sampling. 


For  the  case  of  random  noise,  it  looks  like  the  direct  measurement  method  under 
the  high  rate  sampling  environment.  The  second  method  was  initially  chosen  for  this 
research. 

3.  Method  3.  :  Pre-subtracted  Single  Channel  Measurement  with  Arbitrary 
Sampling  Rate.  (Less  than  Nyquist  Frequency) 

During  the  laboratory  work,  it  was  shown  that  an  unknown  critical  noise  is 
added  to  the  left  sampling  channel  of  the  DPO  which  makes  the  second  method  impos¬ 
sible  to  use.  This  may  be  caused  by  mispositioning  of  the  actual  data  in  the  bulfer.  I  his 
phenomena  and  the  equivalent  model  are  illustrated  in  Figures  5  and  6.  An  alternative 
method  which  avoids  this  problem  is  to  subtract  V  from  A'  using  the  built  in  ".Add  and 
"Invert"  functions  in  the  DPO.  This  results  in  only  one  sequence  of  pre-subtracted  data 
and  is  followed  by  computing  the  estimation  of  the  mean-square  value  of  this  data  se¬ 
quence. 

Referring  to  Figure  6.  assume  that  the  impulse  responses  of  the  sampler  chan¬ 
nels  of  the  DPO  are  unit  impulse  responses,  (  The  effect  of  the  sampler  responses  of  the 
DPO  will  be  discussed  in  the  later  section  of  this  chapter.) 

f(!)  =  ,/(/)  -  o(r)  <2.1 1) 

Then  the  measured  sequences  XP(n)  and  }’,(«)  can  be  denoted  as 

-Vp(n)  =  A’(ii)  +  Dx{n)  4-  .Y»  (2.12.1) 

s  }’(«)  +  Dy(n)  +  Xy{it)  (2. 1 2b) 

where, 

Dx{n)  :  DC  bias  due  to  the  DPO  vertical  adjustment. 

•\»  :  other  noise,  (thermal  noise,  jitter,  quantization  noise  etc.) 
with  similar  representations  for  /),(«)  and  Xy(n)  usig  the  F-channcl. 


A.  Method  1 

N  samples 


AT-l-n 


«x>.(»)=Y  Y  mnk+n) 


TO— 


Delay  Sampling  period 

Dual  channel  measurement  with  sampling  rate  greater  than  Nvquist  freauency. 


B.  Method  2 


N  samples 


i  i 

i - r 


Deiay 


H  T  h- 

Sampling  period 


K- 1 

Rxyin)  =  -jr  Y AUL)  Y(kL  +  n) 
TS  =  LT 


Dual  channel  measurement  with  sampling  r3te  less  than  Nyquist  frequency. 
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C.  Method  3 

N  samples 


-»jnTs}«*- 

Delay 


H  T  re¬ 
sampling  period 


( [  Z{oo)  —  l  Z(ft)  l2]  , 


Z(n)  s  Xp(kL)  -  YJkL  +  n) 


Smgle  channel  measurement  with  sampling  rate  less  than  Nyquist  frequency. 


Figure  4.  Comparison  of  The  Three  Estimation  Techniques  for  The  Crosscorre¬ 
lation 
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Let  a  subtracted  sequence  be  denoted  as 

Z(«)  s  Xp(kL)  -  Yq{kL  +  n)  (2.  i  3) 

where  “ri‘  is  due  to  the  time-shift.  Then  the  mean-square  of  Z(n)  is  described  as 

[  ZUn  r  =  (  X.\kL)  -  yJkL  +  ,i)  |2 

=  [  {  A(kL)  +  Dx(kL)  nx{kL)  {  —  {  1  (k L  +■  n)  -f-  D  a  kL  +//)+■  /jv{kL  +■  n)  J 
=  DC  term  —  2Rxy(u ) 

=  [  Z(oo)  ]*  —  2Rxy{n) 

where  we  have  use umcd  stationarity  of  the  random  process  in  .V  and  V  .  Therefore,  the 
desired  estimation  of  the  crosscoi  relation  function  is 

«.,,(«) -y  [[Z((X)  ]1~  1  Z(l,l  I2  i  ,11J) 

I  he  results  of  the  trial  measurements  reveal  that  this  estimate  has  a  better  and  more 
consistent  value  than  the  result  from  method  2.  This  measurement  technique  has  been 
used  for  this  research. 

C.  PROPERTIES  OF  THE  ESTIMATED  CROSSCORRELATION 

Since  the  crosscorrelation  function  has  an  important  role  in  the  analysis  of  linear 
systems  with  random  inputs,  an  important  practical  problem  is  that  of  dctei  mining  the 
quality  of  the  estimation  of  the  function  for  experimentally  observed  random  processes. 
In  order  to  evaluate  the  quail./  of  this  estimate  it  is  necessary  to  determine  the  mean 
and  the  variance  of  R,v[n).  since  the  estimated  crosscorrelation  is  a  random  variable 
whoso  precise  value  depends  upon  the  particular  sample  function  being  used  and  the 
particular  set  of  samples  taken.  The  mean  value  for  Equation  (2.S)  and  (2.14)  can  be 


AT-! 

£[/U,»)l  -  E[  Y  ^ \kL)  Y(kL  +  n)  ] 

k^  0 
A;- 1 

k  £3 

~jr  K  Rxy{n) 

£| Rk.,(ii)]  =  Rx>(n)  (2.1*i 

rims,  the  expected  value  of  the  estimate  is  the  true  value  of  the  erosscorrelation  function 
and  is  an  unbiased  estimate  of  the  crosscorrelation  function.  The  variance  is  denoted 
as 


=  £|  ).A\v(;m!:1  -2£[  kxy[n)  Rx,{,i)  \  f  £[  {Rx>in)y  \ 

=  £[{axv(«){2i-:  £,>(")  i2 

A -I 

=  £[  {  4*  V AU/.)T(/<£  -h  n) I  -  {RJn)  }2 

/V  ^ 

iisM  0 


A-1  A- 1 


F<f,[  £„(«)  1  =  ~T  Z  Z  ff|{  W)  >V£  +  «) }  (  AT(*£)  Y(kL  +  «) }  J 


Z-<  ^ 

'v  y«o  a=o 
-  (  /<«(»)  }2 


(2.  In) 


Since  j  and  k  arc  different  samples  in  the  ensemble  space,  it  is  reasonable  to  assume  that 
they  are  statistically  independent  random  variables  when  j  *  k.  Hence.  liquation  (Z.lbi 
becomes 

£[!  A VL)Y{iL  +  n) }  \  A  UL.)Y[kL  +  »)}]  =  £|{  X[Kl.)Y[kL  +  /.-) |  .  ./  =  k 

=  A\v("> 

■x  y 

»A 
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,  j*k 


E ({  X{jL)  Y(JL  +  n) }  {  X(kL)  Y(kL  +  «)})={  E\X{kL)  Y<kL  +  n)\  }2 

s  {/?r>.(/d}2 
=  D 

where  A  is  the  crosscorrelation  of X*  and  P  while  B  is  the  square  of  the  crosscorrelation 
for  .V  and  Y.  Using  this  result,  the  variance  leads  to 

l  >[  Rr,tn)  I  =  -jr  (  A.'.  I  +  (  K:  -  K  )  B2 1  -  B1 

l\ 

.1  -  B2 
K 

l  ‘ar[  Rxy(n)  ]  =  —  a2xy(n)  (2. 1 7) 

where  friJ.'O  is  the  true  variance  of  the  product  of  the  random  variables  .V,  and 

An  interesting  phenomena  takes  place  when  t!ie  time  goes  to  infinity.  In  this  cere. 
ii’A’i/t)  and  }'(//)  are  zero  mean  random  variables,  then 

Rxy.(  c«)=.V2F2 
Rx>(00)  =  0 

Therefore, 

r2  '**2 

rtn|  A  Jco)]  =  ^t-  (2.18) 

This  means  that  the  variance  of  the  measured  crosscorrelation  approaches  a  constant 
value  as  the  time  separation  goes  to  infinity.  Therefore,  the  resultant  estimated  cross- 
correlation  could  be  thought  of  as 

RrJn)  =  Rxv(n)  H - ~p — (2.19) 

r\ 

where.  <Hn)  is  the  estimation  error  with  unit  variance.  However,  this  elfc.'-  can  be  in¬ 
duced  by  increasing  the  number  of  sample  points,  since  the  variance  is  norma  ced  In¬ 
dus  number. 


i 


D.  SYSTEM  IDENTIFICATION  USING  ESTIMATED  CROSSCORRELATION 
FUNCTION. 

To  measure  the  input  and  output  noise,  some  kind  of  measuring  device  must  he 
used.  Therefore  it  would  be  reasonable  to  think  that  the  measured  sequence  results  from 
a  convolution  of  the  sequence  to  be  measured  and  the  system  response  of  the  measuring 
device.  This  is  illustrated  in  Figure  7.  Therefore,  the  measured  crosscorrelation  function 
is  Rv  (/)  rather  than  R,y(t). 

R  v,  „<  • )  s  I  -vpt  rh^t  r  +  rtt/r  ( 1 20) 

J  — JO 

where. 


,»oo 


=  Pr)  *  A’(r) 


>Vr)sf  </(.“)  T(t  -  fi)  (in 
=  </(  r)  *  Y(r) 

Substituting  Equations  (2.21)  and  (2.22)  into  Equation  (2.20)  yields 


1*00  j»CO  Q© 

=  p(n)x{T  -  a)  da  \  qiq)y(t  +l  -  n)  df.idv 

J  — oc  J  — O©  J  — OO 

J*  o©  fjz 

-X(T  -  tj)y(  r  +/  -  /O  </r  ,/(/()  d/.uio 

—  >0  J  — .  x?  J  — -V3 

—  f  /’(<?)  f  |  .v(a)j((  /  +  <7  —  //)  +  ,v.)  </j:  </(/.« )  dud  a 

<J  —  j©  J  — y  J  -oc 


R./r  +  «r  —  n)qiii)  dq.da 


12.211 


(2.22) 


) 


Since  the  inner  integral  of  the  right  expression  in  Equation  (2.2')  is  the  convolution  be¬ 
tween  R„(m)  and  </(/i),  the  equation  can  also  be  denoted  as  follows. 


(2.24) 


Ks 


p  OO  pOO 

=  /<<0  q(t  +  o  —  n)Rxy(n)  d\xda 

J  — OO  J  —  OO 

Oyo 

p(a)q((t  -  n)  +  a)  do  Rxy{n)  dn 

-OO 

=  [  Rpq{t  -  P)RXy(P)  dp 

J  — OO 


Therefore,  the  resultant  measured  crosscorrelation  is 


Rx (/)  =  R  [t)  *  Rxv(t) 

<1  r  1 

=  RN(t)  *  /?,(/)  *  h(t) 


•>  ''Si 


=  Rx  AO  *  Kn 

r  <7 


where  is  referred  to  as  the  system  crosscorrclation  in  this  thesis  and  has  the  same 

role  as  the  autocorrelation  in  liquation  (2.?). 

Since  the  measurement  of  the  system  crosscorrelation  is  performed  without  the  tar¬ 
get.  it  should  be  an  even  (unction  if  the  responses  of  the  two  sampling  systems  of  the 
DPO  are  exactly  the  same.  In  fact,  the  responses  of  the  two  channels  do  not  have  exactly 
the  same  characteristics  so  the  test  result  of  the  measured  system  crosscorrelation  has 
an  almost  even  shape,  as  expected.  This  will  be  described  in  chapter  4. 

The  Fourier  transform  of  Equation  (2.25)  results  in 


sv.  (/)  =  s„(/)  sxy(j) 

=  SpqiJ)  sx(j)  //(/) 
=  sXJM)  IHJ) 


(2.26) 


which  satisfies  the  crosscorrclation  theorem. 

To  get  the  impulse  response.  Rr.  (i)  must  be  deconvolved  by  Rt ,  (r).  Iherefoic  at 
least  two  measurements  must  be  performed  for  doing  lime  deconvolution: 


R.xr{»)  =  Rpq(")  *  R.X(n) 


(2.27) 


/?V(I//)  =  f)  *  RJti)  *  hilt) 


I2.2S) 


The  power  spectrutns  of  Equations  (2.27,28)  are 
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(2.29) 


SXrX'(k)  =  Spq{k)Sx(k) 


Sx^(k)  =  Spq(k)  Sx(k)  fl(n) 


(2.30) 


Then,  the  frequency  response  of  the  system  is  extracted  by  dividing  Equation  (2.30)  by 
Equation  (2.29). 


IHk)  = 


Sv ..  (A) 

■V-r 


S.x„xlk) 


(2.31) 


Sn(k)  Sx(k)  !I(k) 


SJk)  SJk) 


However.  this  frequency  division  must  be  compensated  since  the  estimated  cross- 
correlation  function  is  corrupted  by  noise,  as  described  in  Equation  (2.19).  so  that  the 
power  spectral  density  is  biased  by  the  noise  power.  1  he  technique  used  here  lor  noise 
compensation  is  known  as  Riad's  Method.  The  "optimal  estimator"  prescribed  by  Riad's 
Method  [Ref.  10]  is 


//(*)  = 


•V,(*>5W*) 


\Sxx(k)\2  +  klVx 


(2.32) 


where  IV,  represents  the  average  spectral  power  density  of  the  system  crosscorrelation 
function: 


K  Lj 


ik) 


(2.33) 


*=0 


The  parameter  k  is  denoted  as  the  smoothing  parameter  and  establishes  the  lower  limit 
of  the  denominator,  lienee  reducing  the  noise  ell’cct  on  the  computed  target  response 
when  the  input  signal  spectrum  is  small.  The  impulse  response  is  derived  by  performing 
an  inverse  Fourier  transform  on  the  target  frequency  response 
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A(//)  =  F-l{//(A-)} 


(2.3d) 


Another  important  point  to  keep  in  mind  is  that  we  use  an  FFT  algorithm  which 
has  quite  different  properties  from  the  continuous  Fourier  transform.  Ordinarily,  when 
the  FFT  is  used  to  estimate  the  frequency  response,  windowing  is  performed  to  smooth 
the  result.  However,  applying  a  window  interferes  with  the  circular  convolution  prop¬ 
erty.  Therefore,  in  this  research,  we  will  not  apply  windowing  before  performing  the 

* 

FFT  algorithm.  Instead  of  windowing,  the  crosscorrelation  function  i?, ,(«)  will  be  "zero 
padded"  to  avoid  the  wrap-around  problem  which  comes  from  the  circular  convolution 


property. 


III.  BASIC  EXPERIMENTAL  CONFIGURATION 


A.  DESCRIPTION  OF  DPO  AND  ITS  MODIFICATION 

To  meet  the  goals  of  this  thesis  work,  an  ultra  wide  band  noise  source  and  a  dual 
channel  sampling  device  with  high  sampling  rate  and  precise  synchronization  must  be 
used.  A  dual  channel  synchronous  sampling  Digital  Processing  Oscilloscope  (DPO)  and 
a  broadband  noise  signal  generator  were  used  for  this  purpose. 

The  DPO  is  not  an  ideal  sampling  device  since  it  does  not  sample  at  a  rate  exceeding 
the  Nyquist  frequency.  Even  though  the  sampling  rate  of  the  DPO  is  measured  in  MHz. 
it  can  display  a  periodic  repetitive  signal  having  up  to  12  GHz  in  bandwidth  by  rear¬ 
ranging  the  actual  low  rate  sampled  data  using  a  display  buffer.  Hence  the  signal  dis¬ 
played  on  the  oscilloscope  looks  like  it  has  been  sampled  at  over  the  Nyquist  rate.  The 
sampling  technique  of  the  Tektronix  DPO  sampler  S-6  is  illustrated  in  Figuic  f>.  I  he 
Tektronix  S-52  pulse  generator  produces  a  fast  repetitive  pulse  tiain.  A  time  limited 
output  of  a  lumped  network  can  be  displayed  at  an  apparent  sampling  rate  exceeding  the 
Nyquist  rate  by  using  the  S-52  to  drive  the  network.  This  is  shown  in  Figure  9. 

However,  tins  property  could  not  be  extended  to  a  non-repetitive  signal  since  the 
technique  does  not  reconstruct  the  original  shape  of  the  actual  sequence  by  rearranging 
the  sampled  data  after  it  is  taken  by  the  sampling  head.  Consequently,  we  could  not 
utilize  the  virtual  characteristics  of  repetitive  sampling  of  the  DPO  for  this  research.  The 
broad  band  noise  generator  has  about  13  GHz  bandwidth  which  is  far  beyond  the  single 
shot  sampling  rate  of  the  DPO. 

Using  the  measurement  method  2.  as  described  in  Section  II.B,  one  of  the  signals  to 
be  measured  is  manually  shifted  through  the  required  sampling  interval.  The  measure¬ 
ment  process  treats  the  data  as  samples  of  an  ensemble  space  rather  than  a  time  se¬ 
quence  since  our  intention  is  to  compute  the  value  of  the  crosscorrelation  of  the  shilled 
time  difference.  Under  this  condition,  th-'rc  is  no  restriction  of  the  sampling  rate  and  the 
DPO  can  be  used.  There  is.  however,  the  need  of  an  additional  variable  delay  function 
on  one  channel  to  make  the  required  time  shift.  This  function  is  carried  out  using  the 
"Delay"  knob  on  the  channel  0  sampler. 


Figure  8.  The  Sampling  Scheme  of  the  Tectronix  7854  DPO 


A.  A  portion  of  a  rectangular  pulse  produced  by  the  S-52  pulse  generator 
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B.  Modified  pulse  filtered  by  the  DMR005  bandpass  amplifier 
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Figure  9.  A  Pulse  Signal  and  Its  Convolved  Output  by  a  Linear  System 


This  delay  knob  has  a  10  nanosecond  delay  range  which  is  enough  lor  measuring  the 
crosscorrelation  function  for  the  small  targets  that  are  used.  Assuming  that  the  band¬ 
width  of  the  sampled  signal  is  less  than  12  GHz.  the  required  sampling  interval  is 


rs  =  -4-  =  40  psec 

Js 

where,  according  to  the  sampling  theorem, 
/  =  28  =  25  GHz 


on 


0.21 


To  obtain  this  type  of  resolution  on  the  delay,  a  10  turn  potentiometer  was  used  instead 
of  the  built-in  one  turn  potentiometer  on  the  10  nanosecond  delay  knob.  This 
potentiometer  has  a  lice  thousand  division  scale  on  its  face  within  the  full  range  of  its 
delay.  This  modification  is  illustrated  in  Figure  10. 

During  the  laboratory  measurements,  it  was  found  that  the  actual  time  delay  was 
not  exactly  linearly  propottional  to  the  delay  reading.  A  potentiometer  reading  table 
was  then  used  to  mark  the  value  of  the  reading  corresponding  to  each  time  delay.  A 
sample  potentiometer  reading  table  is  shown  in  Figure  II.  Although  the  actual  delay  is 
dependent  upon  the  potentiometer  reading  table  used  with  each  measurement,  a  reading 
dillerence  of  about  3.0  represents  approximately  a  40  picosecond  delay. 

To  create  a  potentiometer  reading  table,  an  additional  measurement  using  the  pulse 
source  was  performed  prior  to  the  noisc-source  measurement.  Alter  creating  several  such 
tables,  it  was  revealed  that  the  potentiometer  reading  has  about  ±  3  picosecond  reading 
error  due  to  misreading  and  the  line-width  of  the  displayed  pulse  when  it  is  performed 
on  20  picosecond  per  division  time  scale.  This  is  illustrated  in  Figure  12.  One  factor 
which  increases  this  line-width  is  the  jitter  noise.  This  reading  error  increases  as  the  time 
scale  increases  and  makes  the  potentiometer  reading  table  next  to  useless.  Therefore,  a 
linear  increment  of  the  potentiometer  by  1.5  is  usually  acceptable  if  the  total  time  win¬ 
dow  to  be  measured  is  over  2<x>  picoseconds.  In  such  a  case  the  measurement  can  be 
performed  using  a  20  picosecond  per  division  time  scale. 

Another  problem  encountered  during  the  laboratory  work  is  a  scaling  dillerence  of 
the  sampling  interval  of  the  two  channels  as  well  as  the  drift  of  the  sampling  frame.  A 
comparison  of(.V—  on  the  2»>  psec.  div.  and  fo  nsec.div.  is  shown  in  I  iguie  13. 

Looking  at  Figure  13,  wc  see  that  the  upper  part  has  noise  with  a  constant  envelope 
while  the  other  demonstrates  a  noise  with  a  consistently  growing  envelope.  Therefore, 


it  would  be  reasonable  to  think  that  the  source  of  the  noise  in  the  upper  figure  comes 
from  the  jitter  noise  of  the  two  channels  due  to  the  random  sampling  time  error  while 
the  other  includes  the  jitter  noise  and  a  scaling  diirerence  of  the  sampling  interval.  One 
possible  source  of  the  scaling  problem  might  be  come  from  machine  accuracy.  So,  if  the 
two  signals  are  aligned  at  one  position  using  the  potentiometer  delay  knob,  so  rhat  the 
subtracted  signal  at  the  point  shows  a  minimum  envelope,  then  the  error  caused  by  dif¬ 
ferent  sampling  rates  grows  with  increasing  the  time  differences  from  the  aligned  point, 
producing  larger  error  in  the  envelope.  Therefore,  it  would  be  better  to  sample  the  data 
on  a  20  picosecond  per  division  (nominal  value)  to  minimize  problems  introduced  by  the 
difference  of  sampling  rates. 

During  the  laboratory  w  ork  it  was  found  that  the  images  of  a  pulse  signal  as  dis¬ 
played  on  each  of  the  two  channels  was  gradually  drifting  from  left  to  right  with  dilleriug 
velocity.  This  drift  problem  might  come  from  the  thermal  inconsistency  of  the  DPO 
since  it  can  be  reduced  by  warming  up  for  about  one  hour.  However,  the  effect  still  re¬ 
mains.  This  problem  requires  completion  of  the  experiment  in  as  short  a  time  as  possi¬ 
ble. 

B.  NOISE  SOURCE  GENERATOR  HARDWARE  AND  IIS  DEVELOI’MEN  T 

A  model  D MR-005  broadband  noise  generator  was  used  lor  this  research.  This 
noise  source  includes  a  noise  generator  and  two  cascaded  broadband  amplifiers,  each  of 
which  has  about  20  dB  gain  with  a  13  GHz  bandwidth.  The  system  output  is  approxi¬ 
mately  400  mV  in  rms  value  and  the  voltage  histogram  resembles  a  Gaussian  drstiib- 
ution.  The  block  diagram,  the  frequency  characteristic,  and  a  sample  recorded  output 
of  the  noise  source  are  illustrated  in  Figures  14.  15,  and  16. 

In  the  laboratory  it  was  found  that  a  frequent  change  of  the  physical  environment 
may  lead  to  electrical  inconsistency  lesuiting  in  output  fluctuations.  A  little  alter  be¬ 
ginning  this  research,  one  of  the  two  amplifiers  used  in  the  noise  source  failed  and  it  was 
sent  to  the  manufacturer  to  fix  the  malfunction.  Alter  the  amplifier  was  returned  to 
NTS.  an  aluminum  support  was  attached  to  the  noise  source  to  avoid  further  problems. 
The  support  also  acts  as  a  heat  sink.  Thereafter,  a  consistent  noise  output  has  been 
produced. 
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Figure  12.  The  Reading  Error  Due  to  the  Line-width  of  Ramp  Signal 


Figure  13.  The  Effect  of  the  Difference  of  Sampling  Interval 
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Figure  16.  A  Recorded  Noise  Signal 


C.  SUMMARY  OF  THE  PROBLEMS  AND  SOLUTIONS  FOR  THE 
CROSSCORRELATION  MEASUREMENT 

The  problems  for  the  crosscorrelation  measurement  and  the  solutions  described  so 
far  are  summarized  as  follows. 

1.  To  produce  a  noise  source  with  consistent  characteristics,  an  aluminum  heat  sink 
support  was  attached  to  the  noise  source. 

2.  A  sampling  rate  lower  than  the  Nyquist  frequency  and  mispositioning  of  the  noise 
added  to  the  channel  0  (output  measurement  channel)  forced  use  of  the  measure¬ 
ment  method  3. 

3.  1  lie  scaling  dilference  of  the  sampling  period  of  the  two  channels  leads  to  sampling 
the  data  on  the  20  picosecond  per  division  time  scale  to  minimize  the  efiect. 

d.  The  drift  of  sampling  frame  requires  a  warm  up  period  of  at  least  one  hour  to  teach 
thermal  steady  state  and  then  to  perform  the  measurement  as  last  as  possible. 

5.  To  perform  the  fastest  measurement,  a  1.5  reading  dilVerence  of  the  potentiometer 
is  linearly  increased  regardless  of  potentiometer  reading  table  since  the  experiment 
reveals  that  a  1.5  reading  difference  of  the  potentiometer  makes  a  20  picosecond 
delay  with  affordable  error  which  can  be  compromised  with  machine  piecisicn  and 
visual  reading  error. 

<>.  Also,  it  is  recommended  to  perform  the  experiment  a  ;  quickly  as  possible  to  reduce 
the  data  transfer  time,  which  is  proportional  to  the  number  o!  the  samples.  (It  t..kc-. 
about  S.5  second  tor  transferring  a  set  of  1024  samples  of  data.) 

Following  these  guidelines,  about  half  an  hour  of  measurement  time  is  needed  for  64 
point  by  2  sets  of  102-1  samples  with  an  additional  one  hour  for  warming  up  the  system, 
which  will  be  the  best  condition  for  the  crosscorrclation  measurement  under  the  current 
laboratory  environment.  However,  the  drift  problem  still  remains  since  at  least  one-half 
hour  is  needed  to  perform  the  experiment. 

Finally,  detailed  tests  revealed  that  the  system  crosscoi  relation  at  the  peak  point  has 
the  value  of  about  twenty  times  the  pcak-to-peak  value  of  the  error  noise.  1  his  is  ''His. 
trated  in  Figure  17.  't  herefore,  if  the  value  of  the  original  crosscori  elation  drops  down 
significantly,  then  the  error  noise  will  dominate  the  measured  signal  and  will  make  the 
crosscorrelation  method  impossible  to  use. 


Figure  17.  The  Amount  of  the  Measured  Error  Noise 


IV.  CALIBRATION  AND  VALIDATION  MEASUREMENT 


A.  INITIAL  TEST  OF  SIMULTANEOUS  CHANNEL  SAMPLING 

The  most  important  and  critical  role  of  the  DPO  in  this  research  is  in  simultaneously 
sampling  a  pair  of  time  functions  which  are  probed  at  the  input  and  the  output  port  of 
a  linear  system.  A  useful  validation  test  for  the  simultaneous  sampling  charactciistic  is 
to  subtract  the  signal  from  one  channel  (.V :  input  sequence  in  this  thesis)  from  the  other 
channel  ( 1’  :  output  sequence).  This  will  be  displayed  on  the  DPO  screen  by  using  the 
built  in  "Add"  and  "Inverse"  functions  of  the  DPO.  Using  the  manual  delay  knob,  the 
image  of  the  subtracted  signal  can  be  minimized  using  the  simultaneous  sampling  posi¬ 
tion.  The  potentiometer  reading  at  this  point  (it  must  be  the  n  =  0  point  of  the  system 
crosscorrelation.  R,^(n))  has  a  value  of  about  200  indicating  that  the  sampling  channel 
t)  i  output  channel)  has  a  time  delay  of  about  2.7  nanoseconds  less  than  that  of  the 
chmncl  1.  This  is  illustrated  in  Figure  IS  with  a  comparison  with  the  original  input  se¬ 
quence. 

Looking  at  the  Figure  IS.  the  left  part  shows  a  noise  with  a  consistent  envelope  of 
small  variance  T  lie  source  of  the  noise  on  the  left  figure  comes  from  the  jittei  noise  of 
the  two  channels  produced  by  random  sampling  time  errors.  Assuming  that  the  envel¬ 
ope  of  the  jitter  noise  is  proportional  to  the  time  derivative  of  the  source  signal  itself,  the 
noise  source  signal  must  have  more  variance  of  its  jitter  noise  vis-a-vis  the  ramp  source 
signal.  This  results  since  the  time  derivative  of  the  ramp  signal  produced  by  DPO  is 
much  smaller  than  that  of  the  noise  source  indicating  that  the  noise  signal  is  more  vul¬ 
nerable  to  the  jitter  noise  than  the  ramp  signal.  This  property  is  verified  in  Figure  FT 

B.  CROSSCORRELATION  OF  THE  DPO  SYSTEM. 

As  described  in  Section  II.D.  the  crosscorrelation  between  input  and  output  is  con¬ 
volved  by  the  crosscori elation  between  the  impulse  response  ol'the  two  sampling  heads. 

I  he  discrete  version  of  eq.(2.2o>  is 


A.  Jitter  noise  of  the  noise  signal 


where. 


k  v,<”)  *  *W")  *  14-) 

Here  /?, ,  (>i)  is  the  estimation  of  the  system  crosscorrelation  and  is  convolved  with  the 

P  4 

crosscorrelation  of  the  sampling  head  and  the  autocorrelation  of  source  noise  itself.  To 
approximately  measure  the  system  crosscorrelation,  we  simply  feed  the  same  input  signal 
to  each  sampling  head,  as  illustrated  in  Figure  20. 

In  the  ideal  case,  each  sampler  system  impulse  response,  pin  and  must  have  the 
same  shape,  and  at  l>est  the  unit  impulse  response.  Although  these  two  impulse  re¬ 
sponses  are  not  exactly  identical,  they  will  in  reality  be  very  close  to  each  other.  Conse¬ 
quently.  the  shape  of  the  system  crosscorrelation  /?  (« I  is  close  to.  but  not  exactly,  an 

even  function,  while  the  autocorrelation  of  the  source  noise  R,(n)  must  be  even. 

The  initial  measurements  were  performed  using  a  2.5  picosecond  sampling  petiud  to 
get  a  more  detailed  time  shape  of  the  system  crosscorrelation.  Subsequent  measurements 
were  performed  using  a  10  picosecond  sampling  interval  to  see  the  emiie  shape  of  the 
system  crosscorrelation.  Figure  21  and  22  show  the  test  results.  As  expected,  the  shape 
of  the  crosscorrelation  is  almost  even,  which  means  the  two  sampling  channels  have 
similar  characteristics. 

C.  DECONVOLUTION  VALIDATION  TEST 

I.  Test  with  Noise  Source 

An  Avantek  model  SA83-2V54  solid-state  microwave  amplifier  was  used  during 
this  validation  test.  This  amplifier  has  an  average  gain  of  42.5  dB  over  the  band  Irom  2 
to  6  GHz  and  provides  a  3  dB  bandwidth  of  about  5  GHz.  The  spectral  characteristic 
of  the  Avantek  amplifier  is  shown  in  Figure  23.  I  his  amplifier  is  linked  to  two  2<»  diJ 
attenuators,  since  it  produces  moderate  power,  to  prevent  any  possible  damage  to  the 
sampling  heads  which  have  an  operating  range  of  5  volts  peak.  Figure  24  illustrates  the 
basic  setup  for  this  test. 

To  extract  the  impulse  response  of  the  Avantek  amplifier  (2-port ».  the  following 
three  steps  were  followed: 

1.  Test  without  the  2-port,  measure  the  system  crosscorrelation. 

2.  Test  with  the  2-port,  measure  the  crosscorrelation  of  the  2-port. 

3.  Perform  deconvolution  computation. 


The  system  crosscorrelation  measurement  was  performed  using  64  time  points 
having  a  40  picosecond  increment.  This  gives  about  a  2.5  nanosecond  window,  which  is 
almost  three  times  longer  than  the  elTective  period  of  the  system  crosscorrelation.  The 
amount  of  the  delay  time  interval  (40  picosecondsi  is  sulliciently  small  in  this  case  since 
this  gives  a  Nyquist  frequency  exceeding  the  bandpass  frequencies  of  the  Avantek  am¬ 
plifier.  A  plot  of  the  measured  system  crosscorrelation  is  shown  in  Figure  25. 

The  bandpass  amplifier  crosscorrelation  measurement  u'as  performed  using  the 
same  time  points  and  delay  interval.  Further,  an  initial  delay  was  given  to  the  channel 
o  in  the  amount  of  400  picoseconds  to  compensate  for  the  effective  delay  through  the 
Avantek  amplifier.  The  resultant  measured  crosscorrelation  is  shown  in  Figure  2<>. 

Finally,  the  computation  of  deconvolution  was  performed  using  Riad  s  optimal 
compensation  technique.  As  described  in  Section  ll.D.  this  deconvolution  method 
serves  to  limit  the  lower  bound  of  the  denominator  and  thus  reduces  the  noise  effect. 
Several  computations  were  performed  with  dilferent  smoothing  factors  and  finally  the 
case  of  •>.!  was  selected  for  the  best  solution.  Alter  the  computation  ol  decomolution. 
the  computed  impulse  response  was  moved  32  points  to  the  right  to  show  the  best 
graphic  interpretation.  Ihe  result  is  illustrated  in  1'iguie  27. 

2.  Test  with  Modified  Pulse  Source 

flic  same  three  measurement  steps  were  performed  using  the  time  domain 
method  described  by  McDaniel  (Ref.  uj.  The  output  ramp  signal  from  the  S-52  pulse 
generator  was  filtered  by  one  of  the  amplifiers  used  with  the  noise  generator  to  ptoduce 
a  roughly  similar  power  spectral  density  at  the  input.  This  source  signal  has  enough 
innate  zero  padding  to  eliminate  the  wrap-around  problem.  The  test  was  fairly  straight 
forward  and  the  results  are  shown  in  Figures  28.  29.  and  30. 

The  comparison  of  the  two  methods  is  given  in  Figure  31.  These  two  results 
have  almost  the  same  shape,  as  is  expected,  aside  from  the  time  scaling  difference.  I  he 
crosscorrelation  method  uses  about  a  40  picosecond  time  interval  while  the  time  domain 
method  samples  with  a  period  of  39  picosecond,  which  is  fixed  by  the  DPO  spaces. 
Another  fact  is  that  the  result  from  the  crosscoi relation  method  has  more  error  noise 
power  than  that  of  the  pulse  method.  Therefore,  the  performance  of  the  crosscoi  relation 
method  using  the  current  DPO  appears  to  produce  tow  er  quality  results  than  that  of  the 
pulse  method. 


Figure  20.  DPO  System  Crosscorrelation  Measurement  Setup 
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Figure  22.  DPO  System  Crosscorrelation  Function  and  Power  Spectral  Density 
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Figure  26.  The  Crosscorrelation  of  Avantek  Amplifier,  Crosscorrelation  Method 
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Figure  27.  The  System  Response  of  Avantek  Amplifier.  Crosscorrelation  Method 
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Figure  28.  The  Modified  Pulse  Input  Signal,  Time  Domain  Method 
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V.  ELECTROMAGNETIC  SCATTERING  MEASUREMENT 


A.  DESCRIPTION  OF  THE  SCATrERING  RANGE 

To  derive  the  characteristics  of  the  noise  source  impulse  scattering  response  meas¬ 
urement,  a  physical  environment  and  equivalent  model  for  the  scattering  range  should 
be  developed.  Since  the  scattering  range  was  already  described  in  detail  by  Mariatequi 
[Ref.  8  ),  and  McDaniel  [Ref.  9],  only  a  brief  description  will  be  given  in  this  section. 

The  general  physical  environment  of  the  scattering  range  is  illustrated  in  f  igure  la. 
As  shown  in  the  illustration,  the  system  can  be  divided  into  5  parts. 

1.  Anechoic  chamber 

2.  Transmitting  and  receiving  antennas 

3.  Broadband  analog  noise  generator 

4.  Dual  channel  sampling  device  (DPO) 

5.  Signal  processor  and  controller  (Tektronix  405 2 A  desktop  computer) 

The  purpose  of  the  anechoic  chamber  is  to  shield  the  incident  and  scattered 
waveforms  from  external  noise  and  interference  effects  as  well  as  multipath  contam¬ 
ination.  The  physical  dimensions  of  the  chamber  arc 

1 .  Longitudinal  length  :  20  ft. 

2.  Lateral  width  :  10  ft. 

3.  Height :  10  ft. 

To  provide  shielding  from  atmospheric  and  man-made  noise,  the  chamber  is 
internally  covered  with  aluminum  panels  which  are  earth  grounded.  A  special  absorber 
material  manufactured  by  Rantcc,  a  division  of  Emerson  Electric,  is  attached  to  the 
aluminum  shield  to  absorb  the  electromagnetic  radiation.  This  absorber  is  ridged  along 
the  length  of  the  chamber  to  guide  the  incident  wave  to  the  back  wall  where  the  most 
of  the  incident  and  guided  electromagnetic  wave  is  absorbed  by  IS  inch  pyramidal  cones. 
Eight-inch  pyramidal  cones  are  attached  to  the  front  wall,  where  the  antennas  are 
mounted.  These  absorb  back-scattered  radiation  and  prevent  rcradiation  of  the 
electromagnetic  wave.  This  material  is  designed  to  provide  back-scattering  attenuation 
of  a  500  MHz  signal  by  about  30  dB  below  that  of  a  Hat  metal  plane.  Reflection  coeffi¬ 
cients  are  increased  to  about  -12  dB  at  1  GHz  for  near  grazing  incidence  on  the  side- 


walls,  lloor  ami  ceiling.  Tliis  leads  to  a  realistic  low  frequency  limitation  of  about  1  GHz 
due  to  the  resultant  multi-path  interference. 

A  styrofoam  column  is  used  as  the  support  for  the  various  targets  in  the  chamber. 
This  material  has  a  relative  permittivity  of  about  1. 1  at  3  GHz  so  that  it  gives  negligible 
elfect  (reflection,  refraction)  to  the  transmitted  wave  but  provides  a  stable  support  for 
the  targets. 

Two  double-ridged  horn  type  antennas  are  used  for  the  transmitting  and  receiving 
antennas.  They  have  a  usable  bandwidth  of  I  to  12.4  GHz  and  have  a  relatively  Hat  gain 
over  the  band. 

The  broadband  analog  noise  generator,  DPO,  and  the  Tektronix  4052A  desk  top 
computer  were  already  described  in  Chapter  III. 

B.  DERIVATION  OF  NOISE  SOURCE  IMPULSE  SCATTERRINC  RESPONSE 
MEASUREMENT  FORMULA 

Since  the  goal  of  this  icsoarch  is  to  validate  the  crosscort elation  measurement  by 
comparing  the  impulse  scattering  response  result  with  that  from  the  time  domain  meas¬ 
urement,  a  simple  sphere  was  chosen  for  the  target. 

An  equivalent  system  diagram  for  the  scattering  range  is  shown  in  f  igure  32.  Here, 
the  impulse  responses  are 

1.  h,{t)  :  Target 

2.  h'{i)  :  Clutter  and  antenna  cross-coupling 

3.  /i,(f)  :  Transmit  antenna 

4.  /i,(/)  :  Receiving  antenna 

and  the  signals  are 

1.  ,r(r)  :  Transmitted  (source)  signal  produced  by  DPO 

2.  ;,(/)  :  Received  (back  scattered)  signal 

3.  .v.( k L )  :  Sampled  transmitted  sequence 

4.  y,(kL  +  n ):  Sampled  received  sequence 

5.  :{n)  :  Stored  sequence  (  A"  -  V) 

6.  R,v(/?|  :  Crosscorrelation  output 

Compared  with  the  discussion  in  Section  IV.C,  the  system  impulse  response  in  the 
anechoic  chamber  is  more  complicated  than  a  simple  2-port  network  response.  The 
impulse  response  of  the  anechoic  chamber  is 


Figure  32.  A  System  Diagram  of  the  Scattering  Range 


(5.1) 


MO  -  Mo  *  (MlO  +  MO)  *  MO 
*  (Mu  *  Mo)  *  (MO  +  hjti\) 

-Mo  *  (Mo  ■+•  Mo) 

This  equation  suggests  the  need  of  at  least  tluee  sets  of  measurements  for  extracting  the 
target  response  MO-  L’sing  our  method,  the  three  crosscorrelations  that  must  be  meas¬ 


ured  are 

/tv.1'0,  -  RpqU)  *  RxyU) 

- 1  Rpq(t)  *  RX{1)  *  MO  *  MO  I  *  Mo  (5.2) 

»  Mo  *  mo 

Rv(n)2  -  0(0  *  (  hpU)  +  hcU )  |  (5.3) 

/M  —  0U)  *  |  M,(  +  h,[i)  |  1 5.4 1 

where, 

0(0  3  /*„(0  *  Rx( 0  *  MO  *  MO  (5.5) 


However,  experiments  revealed  the  need  of  one  more  background  measurement  for  bet¬ 
ter  subtraction  of  the  background  noise.  This  is  because  of  the  use  of  a  relatively  shoit 
time  window  compared  with  the  difference  of  the  propagation  delays  of  dilferent  targets, 
(the  calibration  target,  and  real  target).  Therefore,  instead  of  Equation  (5.2)  to  (5.4), 


we  use  Equation  (5.6)  to  (5.9). 

>  A 

RXpyin)n  =  0(m)  *  M"»2 

(5.6) 

R  («),.  =»  0(«)  *  /!.{«), 

>'7  ' J 

(5.7) 

R.xs,{nh  =  000  *  l  V")  +  M”h  1 

(5.8) 

R-^.  (it)]  =  0(/«)  *  ( /!,(//(  f  /t.OM?  | 

'  5.6 1 

where 

56 


(5.10) 


**■  *1 


4>(n)  *  Rp<j(n )  *  Rx(n )  *  h^n)  *  hr(n ) 


/»,  («) 3  s  hc(n  -  ni)2 


(5.11) 


Here,  m  is  the  time  delay  due  to  the  different  measurement  time  origin  caused  by  the 
physical  shape  of  each  target.  So,  the  best  delay  could  be  determined  by  the  difference 
of  the  transit  time  of  the  leading  edge  of  the  scattered  electromagnetic  wave  for  each 


target. 


•>n 

m  »  Integer  (  } 


In  this  equation.  D  denotes  the  radial  difference  of  the  two  spheres  and  c  represents  the 
free  space  velocity  of  an  electromagnetic  wave.  This  is  illustrated  in  Figure  33  and  an 
example  of  this  effect  is  shown  in  Tigure  34  using  an  I  sq.  ft.  copper  sheet  target  for 
better  reflection. 

The  next  thing  to  do  after  the  lour  measurements  of  the  crosscorrelation  is  the  ex¬ 
traction  of  the  target  response  /i,(/r).  This  is  performed  in  two  steps. 

1.  Clutter  subtraction. 

2.  Crosscorrclation  deconvolution  using  Riad's  method. 


At  first,  subtraction  of  the  clutter  effect  was  performed.  Let  the  subtt acted  result  be 
denoted  as 


3  *Vf(») 2  ~  Kv,(")i2 


=  0(«)  *  hJn) 


(5.12) 


R)  |(«)  = 


(5.13) 


=  </>(«)  *  /tj(rt) 


I  he  follow  ing  is  the  deconvolution  lor  extracting  the  response  of  the  target.  The 
key  point  of  the  deconvolution  is  the  use  of  the  calibration  target  which  has  a  known 
computed  transfer  function.  The  Mic-serics  program  prepared  by  Professor  Morgan 
computes  tire  transfer  function  of  the  sphere  |Ref.  9|. 


HUMBER  OF  SAMPLE  POINT 


NUMBER  OF  SAMPLE  POINT  :  1824 
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Figure  34.  Examples  of  the  Delayed  Signal  Due  to  Different  Paths 


Let  //„,(/)  represent  the  computed  transfer  function  of  the  S  inch  sphere.  For  an  ideal 
case,  the  transfer  function  of  the  target  could  be  computed  as 


m  = 


1 1/) 
Surf 


HpAf) 


(5.141 


where  SM(/)  and  S2l(/)  are  the  Fourier  transforms  of  the  crosscorrelations  R}l(r)  and 
/?,,(r).  Since  we  are  dealing  with  the  estimated  crosscorrelations,  Riad's  method,  as  de¬ 
scribed  in  Section  II.D,  is  applied. 


im  - 


■S3 1 (A)  t | ( k) 

1 5;i(A')j 2  4-  /  ll\x 


uPM) 


(5.15) 


where, 


k=  0 

Finally,  the  smoothed  estimation  of  the  impulse  response  will  be  obtained  from  the  in¬ 
verse  Fourier  transform  of  Equation  (5.15). 

C.  ELECTROMAGNETIC  SCATTERING  CROSSCORRELATION 
MEASUREMENT 

A 

First,  a  12  inch  calibration  target  measurement  of  Rv>(«)2  was  performed.  To  esti¬ 
mate  a  suitable  time  window,  a  pulse  signal  was  first  used  as  the  excitation  in  the 
anechoic  chamber.  The  measured  backscattered  signal  contains  several  constituents:  di¬ 
rect  antenna  coupling;  chamber  clutter  and  target  scattering  which  is  modified  by  the 
antenna  responses  and  multipath.  The  target  pedestal  is  located  approximately  8  feet 
from  the  horn  antennas. 

An  additional  7  feet  of  cable  delay  line  was  linked  between  the  transmitter  power 
splitter  and  the  channel  0  sampler  to  show  the  two  signals  in  the  same  time  window.  At 
first,  a  directional  coupler  was  selected  instead  of  the  power  splitter  since  the  directional 
coupler  transmits  most  of  the  input  signal.  However,  it  was  found  that  the  directional 
coupler  provided  a  significant  attenuation  in  its  high  frequency  band.  The  result  of  the 
alignment  test  is  illustrated  in  Figure  35.  The  leading  edges  of  the  two  signals  arc  aligned 


when  the  potentiometer  reading  indicates  195.  The  time  width  of  the  sphere  response 
was  about  2  nanoseconds. 

Next,  the  noise  source  was  applied  to  measure  the  electromugnatic  scattering  cross¬ 
correlation  measurement.  Since  the  pulse  response  is  limited  to  about  7  GHz  [Fig.  23.], 
a  set  of  64  samples  with  a  sampling  interval  of  40  picoseconds  was  chosen  to  cover  a  2 
nanosecond  time  interval.  This  produces  a  Nyquist  frequency  of  12.7  GHz. 

Unfortunately,  the  measured  value  of  the  crosscorrelation  was  significantly  cor¬ 
rupted  by  the  noise  and  system  errors  so  that  it  was  impossible  to  use.  Since  then,  the 
effort  was  focused  on  finding  the  source  of  the  errors  and  on  the  ways  to  avoid  it. 
However,  the  answer  is  that,  at  this  time,  we  could  not  reduce  the  error  sufficiently  to 
makes  the  estimation  useful.  In  fact,  the  fight  against  the  error  noise  for  the  signal  noise 
source  spans  the  whole  history  of  this  research.  The  situation  in  the  scattering  meas¬ 
urement  offers  signifcantly  less  output  signal  strength  than  was  observed  in  the  initial 
tests  using  lumped  filters.  This  degradation  of  signal  appears  to  have  exceeded  the  limits 
of  fidelity  needed  to  demonstrate  any  viability  of  the  measurement.  In  short,  the  results 
appeared  as  "garbage"  that  could  not  be  reproduced  on  any  two  measurements.  The 
signal  was  essentially  buried  in  the  noise. 

In  addition  to  the  description  about  the  noise  in  Chapter  II,  the  sampled  sequence 
of  the  two  channels  are  modeled  as 

Xp{n)  =  X{n)  *  pin)  +jx{n)  +  Qx{n)  +  dx(n)  +  Nx{n)  +  Dx{n )  (5. 1 7) 

yp(>‘)  -  *'(«)  *  </(«)  +jy{n)  +  Qy(n)  +  dy(n)  +  Sy(n)  +  uy(n)  +  Ny(n)  4-  Dy(n)  (5. 1 8) 

where, 

1.  ,V(m):  Input  signal 

2.  p(n ):  Impulse  response  of  the  input  channel  sampler 

3.  j,(n):  Jitter  noise  of  the  input  channel 

4.  Q,(n):  Quantization  noise  of  the  input  channel 

5.  d,(n):  Drift  (nonstationary)  noise  of  the  input  channel 

6.  Other  noise  (Thermal  noise  etc.) 

7.  [),(")■  DC  bias  due  to  the  DPO  vertical  alignment 

8.  Y(n):  Output  signal 

9.  Impulse  response  of  the  output  channel  sampler 

10.  jy{n):  Jitter  noise  of  the  output  channel 


11.  Q,{n) :  Quantization  noise  of  the  output  channel 

12.  dy{n):  Drift  (nontationarv)  noise  of  the  output  channel 

13.  Sampling  interval  scaling  noise  of  the  output  channel 

14.  uy(n ):  Mispositioning  noise 

15.  ;Y,(/t):  Other  noise  (Thermal  noise  etc.) 

16.  Dy(n):  DC  bias  due  to  the  DPO  vertical  alignment 

It  may  be  surmised  that  the  poor  performance  is  probably  caused  from  the  nonstation¬ 
ary  characteristics  of  the  DPO.  as  summarized  in  Section  1II.C.  This  includes  the  jitter 
noise  and  quantization  noise.  The  estimation  error  for  a  stationary  process  should  be 
reduced  by  increasing  the  number  of  sample  points.  As  shown  in  Figure  36,  the  error 
noise  is  not  linearly  reduced  by  increasing  the  number  of  the  samples.  This  reveals  the 
existence  of  other  sources  of  the  error  such  as  nonstationary  noise  (caused  by  the  drift 
problem)  which  may  give  a  sudden  level  of  noise  that  can  not  be  reduced  by  the  opera¬ 
tor. 

Figure  36  reveals  that  the  electromagnetic  scattering  crosscorrelation  measurement 
for  the  12  inch  sphere  is  significantly  corrupted  by  the  error  noise  although  the  result 
came  from  the  computation  of  8192  samples.  Therefore,  it  might  be  almost  impossible 
to  perform  deconvolution  under  such  a  low  SNR  environment  unless  other  ways  can  be 
found  to  increase  the  SNR  of  the  estimated  value.  One  such  way,  which  may  solve  this 
problem,  is  the  use  of  a  new  generation  of  sampling  devices,  such  as  the  Hewlett  Packerd 
Model  54120  digitizing  oscilloscope.  This  DPO  has  a  sampling  rate  up  to  20  GHz  and 
offers  computer  controlled  time-olfsets  of  its  4  channels.  A  follow-on  investigation  will 
consider  the  application  of  the  IIP-DPO  to  this  research. 
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Figure  36.  The  Fluctuation  of  the  Estimated  Crosscorrelation  Due  to  the  Error 


VI.  CONCLUSIONS 


A.  SUMMARY 

The  main  effort  described  in  this  thesis  is  to  demonstrate  the  viability  of  impulse 
response  measurements  using  a  broadband  random  noise  source.  An  application  of  this 
system  is  radar  target  detection  and  tracking,  where  it  would  offer  the  following  potential 

advantages: 

1.  Low  cost  noise  sources  arc  available. 

2.  Deception  in  radar  applications,  (jammer  or  natural  noise  appearance). 

The  lirst  step  of  the  research  was  to  develop  the  theory  of  noise  source  impulse  re¬ 
sponse  measurements.  An  analog  version  of  the  crosscorrelation  was  derived  followed 
by  development  of  the  three  crosscorrclation  measurement  methods.  T  he  third  meas¬ 
urement  method  (pie-subtracted  single  channel  measurement  with  arbitral y  sampling 
rate)  was  selected  alter  much  effort  in  trying  to  use  method  2.  The  bias  and  variance  test 
of  the  estimated  crosscorrelation  followed.  It  was  shown  that  the  estimation  is  unbiased, 
and  consistent.  Finally,  the  effect  of  the  sampler  on  the  signal  as  well  as  the  application 
of  an  optimal  deconvolution  technique  were  considered. 

The  next  step  was  devoted  to  the  experimental  configuration.  The  characteristics 
of  the  Tektronix  DPO  and  noise  generator  were  investigated.  It  seems  that  the 
Tektronix  DPO  is  not  the  ideal  device  for  this  research  since  it  exhibits  many  problems 
for  these  measurements.  Some  of  these  are  the  drift  problem,  the  relatively  large  amount 
of  jitter  noise  (about  1.  5  of  the  signal  strength),  a  sampling  interval  scaling  problem,  etc. 
An  effort  to  minimize  this  error  noise  was  attempted.  Detailed  tests  were  performed,  and 
it  was  revealed  that  the  system  crosscorrelation  included  errors  in  its  pcak-to-pcak  value 
of  about  1  20  of  its  maximum  signal  power. 

The  calibration  and  validation  test  were  then  performed.  An  Avantck  solid-state 
microwave  amplifier  was  selected  as  the  test  item.  The  crosscorrelation  measurements 
of  the  sampling  system  and  the  Avantck  amplifier  were  pei  formed.  A  smoothed  ftc- 
qucncy  response  of  the  Avantck  amplifier  was  computed  using  the  Riad's  method. 
Then,  the  impulse  response  of  the  Avantck  amplifier  was  obtained  using  the  inverse 
Fourier  transform.  This  was  compared  to  results  obtained  via  a  diiccl  pulse  excitation 


measurement. 


Finally,  an  electromagnetic  scattering  measurement  was  attemped  after  considering 
the  requisite  mathematical  foundation.  Four  measurements  were  needed  in  practice  al¬ 
though  the  theory  requires  three  measurements.  The  scattering  measurement  provided 
insufficient  fidelity  due  to  hardware  limitations  of  the  Tektronix  DPO.  Major  sources 
of  error  in  the  current  system  include  nonstationary  drill,  jitter.and  mispositioning. 

B.  FUTURE  CONSIDERATIONS 

The  main  difficulties  encountered  in  this  research  resulted  from  the  use  of  a  less  than 
ideal  sampling  device  that  could  not  sample  sullicientlv  fast  to  satisfy  the  Nyquist  crite¬ 
rion.  Further,  it  has  been  revealed  that  this  sampling  device  has  many  problems  as  de¬ 
scribed  in  the  previous  section. 

The  significant  payoff  of  this  research  mandates  that  it  be  continued  with  vigorous 
attempts  o  Jleviate  the  hardware  deficiencies  encountered  in  this  initial  efi'ort.  The  use 
of  a  new  generation  digital  processing  oscilloscope  by  Hewlett  Packard  is  a  step  in  the 
right  direction  since  this  device  allows  computer  control  of  the  temporal  ollset  of  each 
channel  with  an  accurate  resolution  of  10  picoseconds.  The  essentially  drilt-lree  per¬ 
formance.  coupled  with  automated  timc-olTscts  in  the  crosscorrclation  data  acquisition, 
will  allow  large  ensemble  averages  to  be  employed  with  much  shorter  measurement  times 
than  have  been  the  case  here. 


APPENDIX  COMPUTER  PROGRAM  LISTINGS 
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